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We present powder and single-crystal neutron diffraction and bulk measurements of the Kagome- 
staircase compound Nia V2O8 (NVO) in fields up to 8.5T applied along the c-direction. (The Kagome 
plane is the a-c plane.) This system contains two types of Ni ions, which we call "spine" and "cross- 
tie". Our neutron measurements can be described with the paramagnetic space group Cmca for 
T < 15K and each observed magnetically ordered phase is characterized by the appropriate irre- 
ducible representation(s). Our zero-field measurements show that at Tph = 9. IK NVO undergoes a 
transition to an incommensurate order which is dominated by a longitudinally-modulated structure 
with the spine spins mainly parallel to the a-axis. Upon further cooling, a transition is induced 
at Thl = 6.3K to an elliptically polarized incommensurate structure with both spine and cross-tie 
moments in the a-b plane. At TLc = 4K the system undergoes a first-order phase transition, below 
which the magnetic structure is a commensurate antiferromagnet with the staggered magnetization 
primarily along the a-axis and a weak ferromagnetic moment along the c-axis. A specific heat 
peak at Tec — 2.3K indicates an additional transition, which we were however not able to relate 
to a change of the magnetic structure. Neutron, specific heat, and magnetization measurements 
produce a comprehensive temperature-field phase diagram. The symmetries of the two incommen- 
surate magnetic phases are consistent with the observation that only one phase has a spontaneous 
ferroelectric polarization. All the observed magnetic structures are explained theoretically using 
a simplified model Hamiltonian, involving competing nearest- and next-nearest-neighbor exchange 
interactions, spin anisotropy, Dzyaloshinskii-Moriya and pseudo-dipolar interactions. 

PACS numbers: 75.25.+Z, 75.10.Jm, 75.40.Gb 
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I. INTRODUCTION 



Quantum spin systems with competing interactions 
can have highly-degenerate ground-state manifolds with 
unusual spin correlations. Small, otherwise insignificant 
perturbations can then become decisively important by 
removing the degeneracy of the low-energy spin fluctua- 
tions, leading to unexpected ground states at low tem- 
peratures. The proximity to quantum phase transitions, 
which separate the various ground states, can lead to new 
types of instabilities which involve charge and lattice de- 
grees of freedom. Examples include exotic superconduc- 
tivity, heavy-fermion conductors and ferroelectricity. 

Frustrated low-spin magnets are ideal model systems 
for the study of competing quantum phases because they 
naturally contain competing interactions in a clearly- 
defined geometry. An important model system is the 
Kagome lattice which consists of corner-sharing triangles 
of spins with equal antiferromagnetic (AF) coupling be- 
tween nearest neighbors. Theoretically it is expected that 
the S = 2 Kagome lattice does not have long-range order 
at zero temperature, but adopts a quantum spin liquid 



ground statei^ The most well known Kagome-related 
magnet is the S = § compound SrCrgGaaOig. However, 
the structure is actually better described as a Kagome 
bilayer with an intervening triangular latticed and it has 
a ground spin-glass like ground state4£ Work on jarosite 
systems showed various types of commensurate long- 
range order stabilized by inter-layer and Dzyaloshinskii- 
Moriya (DM) interactions^! These results indicate that 
Kagome related magnets are highly sensitive to relatively 
weak interactions, and hence arc a likely venue for new 
and exotic ordered states. 

M3V2O8 (NVO) is a system of weakly-coupled spin-1 
staircase Kagome planes containing two types of inequiv- 
alent bonds and magnetic sites. It is thus a variant of the 
highly-frustrated pure Kagome lattice. However, these 
deviations from ideal Kagome geometry introduce sev- 
eral new interactions which determine how the quantum 
degeneracy and frustration are resolved. In addition, it 
has been found 8 that the magnetic ordering generates fer- 
roelectricity and these smaller interactions play a crucial 
role in this phenomenon and can explain the microscopic 
origins of multiferroics. 9 

Magnetic susceptibility and specific heat measure- 



ments reveal that NVO undergoes a series of mag- 
netic phase transitions with temperature and magnetic 
fieldiLii anc [ m pj K ^we reproduce the phase diagram 
as found in Ref. [ll], which we refer to as I. There, our 
neutron diffraction study at zero held showed that NVO 
adopts two different incommensurate phases above 4K, a 
mainly longitudinal incommensurate phase which occurs 
at higher temperature (which we therefore call the high- 
temperature incommensurate, HTI, phase) and a spiral 
incommensurate phase which occurs at lower tempera- 
ture and which we call the low-temperature incommen- 
surate, or LT1, phase. We also found evidence of a com- 
mensurate canted AF phases below 4K. One purpose of 
this paper is to present a comprehensive review of neu- 
tron diffraction data which enables us to characterize the 
HTI, LTI and C phases, but we will not discuss in de- 
tail the C phase. A second purpose is to show that the 
phase diagram can be understood to be the result of com- 
peting nearest-neighbor (nn) and next-nearest-neighbor 
(nnn) exchange interactions, easy-axis anisotropy, and 
anisotropic interactions, both pseudo-dipolar (PD) as 
well as DM interactions. It is particularly important to 
characterize the magnetic structure of NVO in order to 
gain an understanding of the magnetoelectric coupling 8 
which causes the LTI phase to be ferroelectric. 

The experimental details, data analysis, and physical 
interpretation which were only briefly presented in I, are 
fully explained in this paper. We performed a zero-field 
powder and single-crystal neutron diffraction study to de- 
termine the magnetic structures, and we used group the- 
ory to identify the structures that are allowed by symme- 
try for the two observed ordering wave- vectors. Further 
we present the field dependence of the magnetic struc- 
tures by monitoring AF and ferromagnetic (FM) Bragg 
peaks. Magnetic fields up to 8.5T were applied along the 
crystallographic c-direction. We find that application of 
a magnetic field along the c-axis favors the AF C phase 
at the expense of the incommensurate phases. Thus the 
phase diagram obtained by our neutron diffraction ex- 
periments is consistent with that obtained from macro- 
scopic measurements, some of which are presented here 
and provide additional information on the microscopic 
interactions between magnetic ions. 

Almost all our results can be understood on the basis 
of theoretical models that are analyzed in detail in a sep- 
arate paper^ which we refer to as II. To avoid undue 
repetition we will here indicate in qualitative terms how 
these models explain the data and we refer the reader to 
II for quantitative details. 

Briefly this paper is organized as follows. In Sec. II 
we summarize the experimental techniques used in this 
work. In Sec. Ill we give our determination of the crystal 
crystal which confirms earlier work by Sauerbrei et alw* 
Section IV contains the magnetic structure determina- 
tions. Here we give a brief summary of representation 
theory, since this forms the basis for almost all the struc- 
ture determinations. We also present magnetization and 
susceptibility data for magnetic fields applied along the 
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FIG. 1: Phase diagram in H-T when H is applied along 
each of the three crystallographic axes. This phase dia- 
gram is based on specific heat data taken as a function 
of T at constant H fRef llll) . The phases are labeled P 
for the paramagnetic (magnetically disordered) phase, HTI 
for the high-temperature incommensurate phase, LTI for 
the low-temperature incommensurate phase, C for the high- 
temperature canted AF phase, and C for the low-temperature 
canted AF phase. The spin structures of these phases are de- 
scribed below. The transition temperatures (in order of de- 
creasing temperature) are denoted Tph, 2hl, Tlc, and Tec • 
For a field along c, the C and P phases have the same symme- 
try, and therefore there is no sharp transition between them. 



three crystallographic directions. In Sec. VI we give a 
theoretical interpretation of the experimental results. We 
obtain rough estimates of many of the microscopic inter- 
action parameters using the more detailed calculations 
presented in II. Finally, our results are summarized in 
Sec. VI. 



II. EXPERIMENTAL 

Powder samples of NVO were made in a crucible us- 
ing NiO and V2O5 as starting materials^ Single crystals 
were grown from a BaO-V205 fluxiil Neutron diffraction 
experiments were performed using a powder sample with 
a total mass of 10 g and a single crystal with a mass 
of 0.13 g oriented with the (h, k, 0) or (h, k, k) crystallo- 
graphic planes in the horizontal scattering plane of the 
spectrometers. 
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Powder neutron measurements were performed using 
the BT-1 high-resolution powder diffractometer at the 
NIST Center for Neutron Research, employing Cu (311) 
and Ge (311) monochromators to produce a monochro- 
matic neutron beam of wavelength 1.5401A and 2.0775A, 
respectively. Collimators with horizontal divergences of 
15', 20', and 7' were used before and after the monochro- 
mator, and after the sample, respectively. The inten- 
sities were measured in steps of 0.05° in the range of 
scattering angles, 26, from 3° to 168°. Data were col- 
lected at various temperatures from 1.5K to 30K to eluci- 
date the temperature dependence of the crystal structure. 
The program GSASU^ was used to refine the structural 
parameters and the commensurate magnetic structure. 
Additional diffraction data and magnetic order param- 
eters were obtained on the BT7 triple axis spectrome- 
ter to explore the magnetic scattering in more detail. 
For these measurements, a pyrolytic graphite PG(002) 
double monochromator was employed at a wavelength of 
2.47A, with 40' collimation after the sample and no an- 
alyzer. 

The single-crystal neutron scattering measurements 
were performed with the thermal-neutron triple-axis 
spectrometers BT7 and BT9, and with the cold-neutron 
triple-axis spectrometer SPINS. The BT7 experiment 
was performed with 60' collimation after the sample, a 
PG(002) analyzer to reflect 13.408 meV neutrons and 
180' beam divergence after the analyzer. The BT9 
diffraction measurements were performed with an inci- 
dent energy of 30.5 meV and 40'-40' collimation around 
the sample. The H-T magnetic phase diagram was de- 
termined using SPINS with an incident energy of 5 meV, 
a Be filter before the sample and 80'-80' collimation 
around the sample. The SPINS diffraction patterns were 
collected with 80'-80' collimation around the sample, a 
PG(004) monochromator combined with a graphite filter 
in the incident beam and a flat PG(002) analyzer set for 
14.7 meV. The beam divergence between analyzer and 
detector was 240'. 



III. NUCLEAR STRUCTURE 
A. Experimental determination of structure 

The NVO structure refinement from BT1 neutron pow- 
der diffraction data was carried out successfully using the 
previously reported structural parameters^ 3 - as initial val- 
ues. In agreement with previous studies, we found the 
structure to have Cmca symmetry (space group No. 64 
in the International Tables for Crystallography 15 ). No 
structural transition was detected for 1.5K< T <300K. 
The structural parameters and selected interatomic dis- 
tances at two temperatures are given in Table [I] The 
symmetry elements of the Cmca space group of NVO are 
given in Table ITT1 Because V has a low coherent neutron 
scattering cross-section, the atomic positions and tem- 
perature factors of the vanadium ions were fixed to the 
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T = 1.5K 


a (A) 




5.92179(3) 


5.92197(3) 


b(A) 




11.37105(7) 
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Ni 3 
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Ni s -0 2 


2 


2.010(2) 


2.013(2) 


Ni B -0 3 
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2.075(2) 


2.078(1) 
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2.006(2) 


2.006(2) 
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2.083(1) 


2.085(1) 


Ni c -0 3 
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2.0599(7) 


2.0597(7) 


Ni s -Ni c 
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2.9330(6) 


2.9331(6) 


Ni c -Ni c 


2 


2.96089(2) 


2.96098(2) 



TABLE I: Structural parameters and selected interatomic dis- 
tances for NVO, measured using the BT1 spectrometer with 
the Ge (311) monochromator and A=2.0775A. Space group: 
Cmca (No. 64 in Ref. HEl) . Atomic positions (expressed as 
fractions of a, b, and c): Ni a : 8e (notation as in Ref. ITBI) (i, 
y, i); Ni c : 4a (0 ); V: 8f (0 y z ); Oi: 8f (0 y z ); O2: 
8f (0 y z ); O3: 16g (x y z ). Bi = 8n 2 (u 2 }, where u% is the 
displacement of atom i from its equilibrium position and (...) 
indicates a thermal average. Also Rp=^2™ \Il — TV II 
where P a and I\ are the n observed and calculated intensi- 
ties, respectively. Rv, P = \/J2 w i{H — I£) 2 / 2 w i{H) 2 where 
the weight is given by Wi=l/af. The sum of least-squares is 
given by x 2 — Wi (\^o I — Wc\) 2 /{ n — m ), where m is the 
number of fitted variables. 



values obtained by Xray diffraction. 13 

To investigate the effect of magnetic ordering on the 
chemical structure, a series of powder patterns were taken 
below 10K. No significant change for the a- axis, unit 
cell volume, and Ni-0 bond distances were observed, as 
shown in Fig. [21 and Table HJ However, reducing the tem- 
perature below 10K leads to an increasing b lattice pa- 
rameter, while the c lattice parameter decreases without 
a change in space group symmetry. 

The strongest temperature dependence is associated 
with the isotropic mean-square displacements for Ni and 
O, which change significantly with the onset of commen- 
surate order at 4K (Fig. This may indicate a weak 
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Er = 


(x,y,z) 


2 c r = (a, y + 1/2, « + 1/2) 


2 b r = {x,y-\ 


-1/2,2+1/2) 


2 a r = (x,y,2) 
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{x, y,z) 


m ab r = (x,y + 1/2,2 + 1/2) 


vriacX — (x,y 


+ 1/2,2 + 1/2) 


m 6c r = (x,y,z) 



TABLE II: General positions within the primitive unit cell for 
Cmca which describe the symmetry operations of this space 
group. 2 Q is a two-fold rotation (or screw) axis and m a p is a 
mirror (or glide) a — /3 plane. 
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FIG. 2: Lattice parameters a, b, and c as a function of 
temperature, determined from neutron powder diffraction us- 
ing BT1 with a Cu (311) monochromator and a wavelength 
1.5401A. Also shown is the temperature dependence of the 
unit cell volume. 



coupling of the magnetic and the chemical lattice. How- 
ever, we did not observe a crystal distortion with the 
given neutron diffraction resolution. Higher-resolution 
xray or neutron diffraction would be needed to look for 
possible space group symmetry breaking associated with 
the onset of commensurate order. 



B. Structural properties 

Here we note some general features of the structure. 
The structure of NVO consists of Kagome layers of edge- 
sharing NiC>6 octahedra. The layers are separated by 
nonmagnetic V 5+ C>4 tetrahedra. There are two inequiv- 
alent crystallographic sites for Ni, denoted as Ni s and 
Ni c . (We will refer to these as "spine" and "cross-tie" 



sites, respectively.) At 15K, the average Ni 2+ — Ni 2+ 
distance within the layers is d\ = 2.94A, while the inter- 
layer distance is d,2 — 5.69A. Based on the relatively 
large inter-layer to intra-layer ratio, difd\ = 1.9, a strong 
two-dimensional magnetic character may be expected in 
this compound, with the magnetism dominated by intra- 
layer Ni 2+ — Ni 2+ exchange interactions. Unlike pre- 
viously studied Kagome lattice-based materials, which 
have planar magnetic layers, the Ni-0 layers in NVO are 
buckled, resulting in the "Kagome-staircase" structure. 
The symmetry of the superexchange interaction medi- 
ated by O ions shows that there are two inequivalent su- 
perexchange paths between neighboring Ni 2+ ions within 
Kagome planes. In particular, there is a superexchange 
path between Ni s positions along the crystallographic a- 
direction which is different from nn interactions between 
Ni ions on neighboring Ni s and Ni c positions. 
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FIG. 3: Isotropic mean-square displacements of Ni and O as 
a function of temperature, obtained from the same powder 
spectra mentioned in Fig.0 

The point group is Abelian and its three generators 
are spatial inversion (X), and the two- fold rotations R a 
and Rb about the a and b axes respectively. The space 
group operations are X and R a referred to the origin at 
a cross-tie site, such as site cl in Fig. 0] and Rb about 
an axis passing through a spine site such as site si in 
Fig. 0] The operations give rise to mirror planes perpen- 
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dicular to the a-axis passing through a cross-tie site and 
a two-fold screw axis about a z-axis passing through a 
spine site such as site si with a translation of c/2 along 
c. There is a glide plane perpendicular to the c-axis pass- 
ing through a chain of spine sites with a translation a/2 
along a. The result of these symmetry operations is that 
all the Ni spine sites are related by symmetry and all 
the Ni cross-tie sites are similarly related by symmetry. 
Our convention for the coordinate axes is as follows: the 
spines he along the (100), or a-axis (denoted sometimes 
the IE-axis), the basal plane includes this axis and the 
(001), or c-axis (or sometimes the z axis), and the axis 
perpendicular to this plane is denoted (010) or the b-axis 
(or sometimes the y-axis) . The two types of Ni sites are 
shown in Fig. 0]and their coordinates are given in Table 

eh 



t A 




H 



FIG. 4: Ni sites in the primitive unit cell, whose vertices are 
A, B, . . .H and is defined by the vectors vi = (o/2)o+ (6/2)6, 
V2 = (a/2)a — (6/2)6, and V3 = cc. The "cross-tie" sites (on- 
line=blue) cl and c2 lie in a plane with 6 = 0. The "spine" 
sites (on-line=red) are labelled si, s2, s3, and s4 and they 
may be visualized as forming chains parallel to the a-axis. 
These chains are in the buckled plane with 6 = ±5, where 
S = 0.136 as is indicated. Cross-tie sites in adjacent planes 
(displaced by (±6/2)6) are indicated by open circles. Spine 
sites in adjacent planes are located directly above (or below) 
the sites in the plane shown. 



The spine and cross-tie sites have different local sym- 
metry and will be seen to have very different magnetic 
properties. In the presence of AF ordering of the spine 
sites, the cross-tie sites are frustrated if one assumes 
isotropic AF Heisenberg interactions; the sum of the mo- 
ments of the four spine neighbors of each cross-tie spin 
vanishes. In this regard, this system is reminiscent of 
Sr 2 Cu304Cl2i& and of various "ladder" systems which 
have recently been studied>i£ 

The above structure has several implications for the 
magnetic interactions. As explained in Ref. Hoi the lead- 
ing nn Ni-Ni magnetic coupling arises via superexchange 



= (0.25,-0.13,0.25) 
r 2 s = (0.25,0.13,0.75) 
= (0.75,0.13,0.75) 
= (0.75,-0.13,0.25) 
rl = (0,0,0) 
vl = (0.5,0,0.5) 



TABLE III: Positions of Ni 2+ carrying S=l within the prim- 
itive unit cell illustrated in Fig. 0] Each component is ex- 
pressed as a fraction of the respective lattice constant, so 
= 0.25a — 0.13b0.25c. Lattice positions r™ are spine sites 
and r™ are the cross-tie sites. NVO orders in space group 
Cmca, so there are six more atoms in the orthorhombic unit 
cell which are obtained by a translation of lattice Vj through 
r 6 by (0.5,0.5,0). 

interactions, mediated by two Ni-O-Ni bonds. For a pair 
of nn spine spins, the angles of these two bonds are 90.4° 
and 95.0°. For a spine-cross-tie pair, these angles are 
90.3° and 91.5°. For the similar case of Cu-O-Cu bonds, 
it has been shown that when these angles are close to 90° 
then the resulting exchange energy is small ^ and may 
even change its sign from FM to AF (as the angle de- 
creases from 90°). Since both Ni and Cu involve d-holcs 
in the high e g states (within the oxygen octahedron sur- 
rounding Ni or Cu) , we expect similar results to apply for 
the Ni case. Accordingly, we do not necessarily expect 
that nn interactions dominate second neighbor interac- 
tions. Similar calculations for the related nnn coupling, 
via Cu-O-O-Cu, gave AF interactions^ Similar Ni-O-O- 
Ni interactions could compete with the nn interactions, 
and give rise to incommensurate structures, as explained 
below. 

As a result of the crystal symmetries, there is a limited 
number of independent nn interaction matrices. If we 
write the interaction between spine spins S(i) and S(j) 
as 

Ha = ^2M aP {i,j)S a (i)S fi {j), (!) 

a/3 

where a — a, b, c is a component label, then once we 
have specified M a p(sl, s4) (in the notation of Fig. 0J, 
we can express the interaction matrices for all other nn 
pairs of spine sites in terms of M Q . | g(sl, s4). Similarly, 
we only need to specify a single interaction matrix, e. g. 
M a p(sl, c2), for nn pairs of spine and cross-tie sites or for 
interactions between nnn in the same spine. Symmetry 
also places some restrictions on the form of M a p(sl, s4). 
These results are obtained in II. 



IV. MAGNETIC STRUCTURE 

A. Representation theory 

Since group theoretical concepts are central to our 
analysis, we shall here summarize the results which we 
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will invoke. Additional details are available in II and in 
Appendix A. In general, Landau theory indicates that 
the free energy F in the disordered phase is dominantly 
a quadratic form in the spin amplitudes S a (r) at site r. 
Thus 

F = \ E «a/3(r,r / )5 Q (r)S /3 (r') . (2) 

r,r' 

In view of translational invariance this free energy can 
be written in terms of Fourier amplitudes as 

F= \ £ »«rJJr'(Q)Sa(Q,T)S / ,(-Qy) ) (3) 
Q,t,t',q,/3 

where Q is the wavevector— and r labels sites within the 
unit cell. Here the Fourier amplitudes are given by 

S a (Q,r) = A^]>> Q (R + r r ) e iQ - (R+rT) , (4) 

R 

where r T is the position of the rth site within the unit cell 
and the sum over R is over translation vectors for a sys- 
tem of N uc unit cells. As our data indicate, the ordering 
transition is a continuous one which is signaled by one 
of the eigenvalues of the quadratic free energy passing 
through zero as the temperature is reduced. This condi- 
tion will be satisfied by some wavevector, or more pre- 
cisely, by the star of some wavevector q, which is in the 
first Brillouin zone of the primitive lattice. (This is usu- 
ally called "wavevector selection.") The critical eigenvec- 
tor, i. e. the eigenvector associated with this instability, 
will indicate the pattern of spin components within the 
unit cell which forms the ordered phase. In view of the 
symmetry of the paramagnetic crystal, which dictates the 
form of Eq. © , we see that this eigenvector has to trans- 
form according to one of the representations of the sym- 
metry group which leaves the wavevector q invariant. 20 
(This group is usually called "the little group.") The as- 
sumption that the instability (toward the condensation 
of long range order) involves only a single representation 
is based on the assertion that there can be no accidental 
degeneracy (which would correspond to a higher order 
multicritical point). Here we are interested in the rep- 
resentations of two types of wavevectors, namely zero 
wavevector (in which antiferromagnetism arises because 
of AF interactions within the unit cell) and an incom- 
mensurate wavevector (q, 0, 0) at some nonspecial point 
on the x-axis. Because of the magnetic structure of the 
unit cell one has to be careful in relating the selected 
wavevector (q, 0, 0) to the physically relevant quantity, 
namely the Fourier coefficient of the magnetic moment. 
The representations for the wavevector (q, 0, 0) are de- 
scribed in an Appendix. Within a given representation T 

(T) (T) (V) 

one has the parameters m sa , m sb an d m-sc which com- 
pletely fix the a, b, and c components, respectively, of 

the Fourier amplitudes S a (Q,sn) of all the spine spins 

(T) (r) 

within the unit cell and the parameters m ca , m cb , and 



nice which similarly completely fix the Fourier ampli- 
tudes S a (Q, cn) of all the cross-tie spins within the unit 
cell. For some representations some of these amplitudes 
may not be allowed to be nonzero. For an incommen- 
surate wavevector these parameters are complex-valued, 
although, of course, the resulting spin components must 
be real, because we should invoke both T associated with 
q and T* associated with q. 

Happily, as the above discussion implicitly assumes, 
the situation is quite simple in that for many systems, 
such as NVO, all the representations are one dimensional. 
What that means is that under any group operation, the 
allowed eigenvectors transform either into themselves or 
into a phase factor of unit magnitude times themselves. 
To summarize the results of Appendix A, if O p is an op- 
eration that leaves the incommensurate nonzero wavevec- 
tor invariant, then we may write 

P mf2 = &CT)m£> , (5) 

where g assumes the values s for spine and c for cross- 
tie, a = a, 6, or c, and £p(T) is the character for the 
symmetry operation O p in the representation T. For zero 
wavevector (relevant for the AF phases) these characters 
are given in Table IVl and for the incommensurate phases 
they are given in Table IVTT1 

The discussion up to now took account only of those 
operations which leave the wavevector invariant. How- 
ever, the free energy must be invariant under all the sym- 
metry operations of the paramagnetic phasepSi* 2 ^ Thus 
far the symmetry properties we have discussed apply to 
any crystal whose paramagnetic space group is Cmca. 
Now we discuss the consequences of restricting the mag- 
netic moments to the spine and cross-tie sites which have 
higher site symmetry than an arbitrary lattice site and 
in particular we will study the incommensurate phases. 
We consider the effect of spatial inversion on the spin 
wavefunctions. To illustrate the concepts we assume a 
wavevector q (in units of 2n /a) along the a-axis and con- 
sider a spin configuration which transforms according to 
T4 and which therefore has the components ipi given in 
Table l\Tfl Since the magnetic moment is an axial vec- 
tor, spatial inversion X takes the moment into itself but 
moves it to the spatially inverted lattice site. Let us con- 
sider the spin state at spine site #3 in the unit cell at 
R = (X, Y, Z) . Before applying spatial inversion the spin 
vector at that site is 

S 3 (R;r 4 ) = [m sa k-m sb b + m sc c}e 2 ^ x+x ^/ a 
+ [m sa k - m sb b + m sc c]*e- 2 ^ x+x ^' a , (6) 

where x s % is the ^-coordinate of spine site s3 within the 
unit cell and the representation label is implicit. After 
inversion (indicated by a prime) the spin at this site will 
be that which before inversion was at — R — x S 3a, which 
is a site of sublattice #1. Thus 

S 3 (R; T 4 )' = [m sa k + m sb h + m sc c]e- 2m ^ x+x ^' a 
+[m sa a + m sb b + m sc c]*e 2 ^ x+x ^' a . (7) 
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By comparing Eqs. I© and J7J) we see that for a 
a = c we have that 

1[m sa ] = [m sa \* , 

and for the ^-component we have 



a or Taking account of Eq. (|13|l this implies that 



Irrist, = —m 



si, 



(8) 



(9) 



One can check that the spin components of the other 
spine sites transform this same way. Furthermore, this 
type of analysis indicates that all the coordinates of Table 
IVIII for the cross-tie sites obey Eq. (JSJ). 

Accordingly, for this irreducible representation (irrep) 
we now introduce symmetry-adapted coordinates rh ga 
which obey Eq. ©. We write 



???.„ 



(10) 



except for g — s and a = b. For this case, to trans- 
form coordinates so that Eq. © is transformed into the 
desired form of Eq. JSJl , we write 



m sb = im sb 



(11) 



Ba 
'J i 

g'g 



aB 

V 



Ba* 
' J 9'9 



(15) 



In other words, for the present symmetry, the coeffi- 
cients in the quadratic free energy (expressed in terms 
of symmetry-adapted coordinates) are all real-valued! 
What this means is that the eigenvector of the quadratic 
form can be expressed as a single overall complex phase 
factor times a vector with real- valued components. This 
condition ensures that all the amplitudes which make up 
the spin eigenvector have the same phase. 

The above discussion incorporates several implicit as- 
sumptions. For instance, it is assumed that the magnetic 
moment is truly localized on the high symmetry Ni sites, 
whereas in reality this moment is spread out around these 
sites. Also, in principle even the "nonmagnetic" oxy- 
gen atoms will have a small induced magnetic moment. 
In addition one can consider the effects of quartic and 
higher-order terms in the Landau free energy as well as 
fluctuations not included by mean field theory. These 
corrections are analyzed in II. 



For the other irreps the analysis is similar, but the com- 
ponents which have to transform as in Eq. (|11|) may be 
different. The complex- valued symmetry adapted coordi- 
nates which transform according to each of the irreps are 
collected in Table I Villi and all of these are constructed 
to obey Eq. ©. 

When we assume the condensation of a single irrep, T, 
the Landau expansion in terms of the above symmetry 
adapted coordinates is of the form 

F = \ E »%K*™*1> > ( 12 ) 

gg'aB 

where fh ga is shorthand for fh T ga and the reality of F 
requires that 



Because these coordinates transform according to a sin- 
gle one-dimensional irrep we know that this expression 
for the free energy is indeed invariant under the oper- 
ations of the little group. But the free energy is also 
invariant under spatial inversion. This additional invari- 
ance provides addition information. This situation has 
been reviewed recently by J. Schweizer— , but the ap- 
proach we use below may be simpler in the present case. 
Here, because of the special transformation property of 
Eq. © we have that 

F = XF= \ E v gg'( Irfl 9^y ir ng'P 

gg'aB 

= \ E v%(™ g *)(m gl0 y . (14) 

gg'aB 



B. Domains 

We now discuss the effect of domains for the case of 
zero applied field. In this case if we assume a single spin 
eigenf unction, ^i, then, because we condense order out 
of the paramagnetic phase, we expect to have all eight 
states of the type O^i, where Oi is one of the eight sym- 
metry elements of the paramagnetic space group. The 
simplest way to take account of these possibly different 
structures is to associate with a given scattering vector 
Q, the intensity averaged over the set of eight wavevec- 
tors OiQ. These eight states will not, in general, be 
distinct, but this is a simple automatic way to take do- 
mains into account. Indeed, for the HTI phase, these 
eight states will generate four times the state *i and 
four times the state — Since these two configurations 
both give the same neutron scattering signal, there is no 
effect due to domains, as is the case for two sublattice 
antiferromagnets. The result is less trivial when we have 
simultaneous appearance of two different representations 
as in the LTI phase, where we have T4 and T±, where the 
Fourier components can either be added or subtracted 
from each other. Within the accuracy of the experiment 
we could not distinguish whether or not both such do- 
mains occur simultaneously in NVO. 

This procedure can be extended to nonzero magnetic 
field H. If the field is large enough (or has been obtained 
by reducing the field from a large initial value), one would 
assume that all domains obtained by applying the oper- 
ations Oi of the subgroup of the space group which leaves 
H invariant occur with equal probability. Then one as- 
sociates with a given scattering vector Q the intensity 
averaged over the set of wavevectors o^Q. 
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C. Magnetic order at zero field 
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FIG. 5: Low-angle portions of the BT1 neutron powder 
diffraction pattern collected at 15K, 5K, and 1.6K. (a) Nu- 
clear structure fitting, (b) Nuclear structure fitting only, in- 
commensurate magnetic peaks observed, (c) Both nuclear 
and magnetic structures were included in the fit. The differ- 
ences between observed and calculated intensities are shown 
at the bottom of each figure. The vertical lines indicate the 
positions of the possible Bragg peaks. 



Fig. shows the low-angle neutron powder diffraction 
pattern measured at 1.6K, 5K and 15K. The appear- 
ance of new Bragg peaks upon cooling indicates that the 
compound undergoes transitions to magnetic order below 
9. IK. The Bragg peaks below 4K can be indexed with 
ordering wave-vectors that are commensurate, whereas 
no such identification is possible for higher tempera- 
tures, suggesting incommensurate magnetic structures at 
higher temperatures. Fig. [fj] shows the temperature de- 
pendence of the intensity of an incommensurate magnetic 
peak near 26 = 21.49°, and scattering associated with 
the commensurate order observed at 20 = 45.05°. This 
shows that the incommensurate phase exists in a finite 
temperature window between 4 and 9. IK, and that the 
low T state has commensurate magnetic order. 
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FIG. 6: Top: Intensity at the scattering angle 2Q, close to 
the AF peaks (1,1,2) and (1,3,0). The intensity above 4K 
is related to (0.73,3,1) magnetic reflection associated with 
the incommensurate magnetic structures. Bottom: Intensity 
at the scattering angle 20 for the incommensurate magnetic 
peak (0.73, 1,0) as a function of temperature. 



The magnetic order was further investigated with ex- 
periments on a single-crystal in which scattering was 
monitored only in the (h, k, 0) and (h, k, k) planes. Fig.0 
shows the elastic neutron scattering at three differ- 
ent temperatures for wavevectors of the form (Q x , 1,0). 
Upon entering the magnetic phase, a Bragg peak is 
formed with a maximum intensity for Q x = Q® w 1.27 r. 
1. u.'s. This result indicates weight in the Fourier trans- 
form of the spin 5(Q) at a wavevector {Q%, 1, 0) which is 
outside the first Brillouin zone of the primitive unit cell. 
We deduce that, Bragg scattering is allowed at wavevec- 
tors 



G + q= G± (g,0,0) 



(16) 



where q 0.27 and G is a reciprocal lattice vector of the 
form G = (l+m,l — m,n), where I, m, and n are integers. 
The ordering wave-vector in the HTI and LTI phases is 
thus v = (q, 0,0). The peak in Fig. 7 is in the I = 
l,m = n = zone. In this formulation the wavevector 
q indicates that the spin function varies as exp(iq • R) 
as the position is displaced through a translation vector 
R of the lattice, as defined in the caption to Fig. 0] 
Thus, for a given irrep, the actual spin amplitudes are 
determined by the value of q and the value of the spin 
coordinates within the unit cell as given in Table IVTT1 for 

Eng. 

Note that this wavevector q does NOT give the phase 
factor introduced upon moving from one spine site to 
its nearest neighbor. In view of the intra-cell structure 



9 



1400 




1.22 1.24 1.26 1.28 1.30 1.32 



FIG. 7: Neutron diffraction intensity measured as a function 
of Q x for wavevector Q=(Q I ,1,0). Here the peak position 
determines the value of 1 + q. 



(given in Table IVTT|l . different components of spin will 
involve different phase factors. If we consider the a- 
component of spin in irrep T 4 and let the translation 
vector R be (X, Y, Z) , then we see that 

S.l, a {X) = e M^/a)[X+(a/i)] m a ; 

S si , a (X) = e «?(2"/a)[X+(3a/4)] ( _ m a) > ^ 

where si and s4 are site labels as in Fig. |2| 
Thus S a 4 ia (X)/S sl>a (X) = -e iqw . Similarly one finds 
S s i, a {X + a)/S S 4,a(X) — —e lqv . Our conclusion is that 
translation along a spine by a/2 introduces a phase fac- 
tor e^ 9+1 - )7r . So, the wavevectors for the a-component 
along a single spine are q± = 1 ± q. As we shall see later, 
other wave-vectors are needed to reproduce the position 
dependence of other spin components within this repre- 
sentation. 

The incommensuration q is weakly temperature de- 
pendent, as shown in Fig. [SJj, indicating competing in- 
teractions in the spin lattice. The temperature depen- 
dence of the integrated intensity of the incommensurate 
Bragg peaks indicates the onset of magnetic order at 
Tph = 9. IK and further a second order transition at 
about Thl = 6.3K. These transition temperatures are 
consistent with specific heat measurements, which show 
sharp peaks at these temperatures, and with magnetiza- 
tion data. The existence of two different incommensurate 
magnetic phases is a further indication of competing in- 
teractions in NVO. 

Fig. [S] shows that the incommensurate Bragg peaks 
abruptly lose most of their intensity at Tlc = 4K, be- 
low which temperature a commensurate magnetic order 
becomes dominant. Commensurate Bragg peaks were ob- 
served at the (h, k, 0) for h + k = even, so the commensu- 
rate structure is associated with an ordering wave- vector 
v = (0,0,0). The magnetic unit cell in the C phase is 
thus identical to the chemical unit cell. 

In the zero-field cooled sample, we observed a weak in- 
commensurate peak which is not present in the 8T field 
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FIG. 8: Temperature dependence of magnetic Bragg peaks in 
zero field. The low-temperature field was reached by zero-field 
cooling. Integrated intensities were obtained by integrating 
diffraction intensities measured as a function of h wave- vector 
transfer, (a) Temperature dependence of /i-integrated inten- 
sity of the AF (1, 1, 0) reflection, (b) Temperature dependence 
of /i-integrated intensity of the incommensurate (1 — q, 1,0) 
reflection, (c) Temperature dependence of the wave-vector 
(1 ± q, 1, 0) of the incommensurate magnetic reflections. (The 
data points below T = 4K represent nonequilibrium domains 
of the LTI phase embedded in the C phase. Such peaks are 
only present after cooling through the LTI phase.) (d) Tem- 
perature dependence of /i-integrated intensity of the incom- 
mensurate (1 + q, 1, 0) reflection. 



cooled sample. This is evidence that the incommensurate 
phase below Tlc is metastable - possibly a reflection of 
how close the commensurate and incommensurate mag- 
netic order lie in energy. However, Fig. [5] shows that the 
ground state of NVO can be annealed through the ap- 
plication of a field along the c-axis which makes the in- 
commensurate Bragg peak vanish. The incommensurate 
Bragg peaks do not reappear when lowering the field, but 
the intensity is instead transferred to the commensurate 
peak that grows in strength. 



D. Field dependence of Bragg reflections 

The field dependence of the magnetic Bragg reflections 
was investigated only with fields along the c-axis. Fig.llUI 
shows the magnetic reflections at 2T as a function of 
temperature. Upon heating, the (1, 1, 0) reflection disap- 
pears in a first-order transition as the incommensurate 
(1 ± g, 1, 0) reflections appear. The commensurate phase 
survives to higher temperatures than at zero field, and 
the ordered moment increases with field, both indications 
that the magnetic field stabilizes the commensurate or- 
der. The LTI magnetic structure occupies a relatively 
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narrow temperature range while the temperature bound- 
aries for the HTI magnetic structure are nearly indepen- 
dent of field. 

The LTI magnetic structure is further suppressed with 
increasing magnetic field along the c-axis. At 5T, as 
shown in Fig. ^JJ the LTI structure does not occur, and 
as the temperature is increased, the commensurate struc- 
ture gives way directly to the HTI magnetic structure. 
The phase transition between the paramagnetic and HTI 
phase occurs at a temperature somewhat below its zero- 
field critical temperature Tph = 9K. It is for this field di- 
rection that the phase boundaries depend most strongly 
on the field. The magnetic phase boundaries obtained 
with these neutron measurements are consistent with the 
phase diagram obtained through specific heat measure- 
ment with the field along the c-axis. 

A field along the c axis in the HTI phase leads to sup- 
pression of the incommensurate Bragg peak at (1 — 1, 0) 
and an increase in the intensity of the (1, 1, 0) reflection, 
as shown in Fig. ^Jfor T = 8.4K. The temperature de- 
pendence of the intensity of the incommensurate Bragg 
peaks suggests that the HTI phase disappears in a contin- 
uous phase transition at a critical field H c , giving way to 
a commensurate field driven AF phase at higher fields. 
This contrasts with the first order nature of the phase 
boundary between the LTI and AF phases. As shown in 
Fig. 1121 both the incommensurate wave vector, q, and 
the integrated intensity of the (1,1,0) AF Bragg peak in- 
crease progressively more rapidly as the field increases. 



E. Phase diagram 

The zero- field phase boundaries at Tph, Thl and Tlc 
deduced from the diffraction experiments are consistent 
with those observed with specific heat measurements. In 
contrast, the intensity of the (1, 1,0) does not show any 
anomaly at Tec, to a level of 0.5%. This suggests that 
the CC phase transition leaves the magnetic structure of 
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FIG. 9: Integrated intensity of the (1, 1, 0) and the (0.72, 1, 0) 
reflections for increasing and decreasing field after the sample 
was cooled in zero field. 



C-phase described by ordering wave-vector v = (0, 0, 0) 
unaltered. 

Figures EH and ^] already show that the C-LTI and 
the C-HTI transitions are first order at small H . Simi- 
lar results were obtained when we varied H at fixed T. 
Specifically, Fig. [21 shows the jumps in the staggered 
moment of the C phase as one moves from the LTI or 
from the HTI phases into the C phase, for T < 8K. In 
contrast, at higher temperatures and fields the transition 
from the HTI to the P phase is continuous, as can be seen 
from Fig. ^| In fact, at finite field along c one cannot 
really distinguish between the P and the C phases, and 
we already saw in Figs. 1101 and HTI that the HTI-P tran- 
sition is continuous for H < 5T. We thus conclude that 
the HTI-C transition changes from being continuous to 
being first order somewhere around the top of the bound- 
ary of the HTI phase in Fig. ^ Such tricritical points 
are abundant in anisotropic antiferromagnets subject to 
magnetic fields. 

The change in slope of the (110) intensity versus T 
curve at the HTI-P transition (see Fig. I12|) indicates 
the coupling between the commensurate and incom- 
mensurate order parameters, which will be discussed in 
Sec.lVDl 

The observations which are relevant for the C phase are 
(a) the specific heat data at 9T which show no anomaly 
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FIG. 10: Temperature dependence of magnetic Bragg peaks 
for H = 2T, applied along the crystallographic c-axis. The 
low-temperature state was reached by cooling the sample in 
a field of 8T. Integrated intensities were obtained by inte- 
grating the intensity of the magnetic Bragg reflection when 
opbserved by a scan over the /i-component of the wave- vector 
transfer, (a) Temperature dependence of integrated intensity 
of the AF (1, 1, 0) reflection, (b) Temperature dependence of 
integrated intensity of the incommensurate (1 — q, 1, 0) reflec- 
tion, (c) Temperature dependence of the incommensuration 
q. (d) Temperature dependence of integrated intensity of the 
incommensurate (1 + q, 1, 0) reflection. 
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and (b) the lack of any detectable anomaly in the In- 
dependence of the(llO) peak intensity (see Fig. IT3|l . As 
will be shown in Sec. IV C II the P phase in a field allows 
the development of order characteristic of only the irrep 
Ti in which there is a uniform moment along c and a 
staggered moment along a. This suggests that the mag- 
netic structure in the C and P phase (in a field) exhibit 
the same symmetry, and that the symmetry at high field 
at low and high temperature is the same as that of the 
C phase at zero field because it can be reached without 
crossing a second order phase transition. 



F. Magnetic Structures 

Even more detailed information about the spin inter- 
actions can be obtained by determining the symmetry of 
the ordered magnetic structures. In this paper, we will 
focus our attention on the HTI, LTI, C and P phases, 
and we will leave a discussion of the C phase for a later 
publication. 



where the number in parenthesis is the uncertainty in 
the last digit quoted. The quality of the fit is given by 
,Rp=0.15 and x 2 =12 (see Table H] for the definition of 
these quantities). The spin arrangement at T — 7K is 
illustrated in Fig. 115b . It is predominantly a modulated 
structure on the spine sites with moments that are paral- 
lel to the a direction. The moments on the cross-tie sites 
either vanish or are very small. [If the cross-tie moments 
in the b direction are nonzero, they are out of phase with 
the spine moments because of the phase factor introduced 
by Eq. (fTTf) ] . 

To determine the HTI magnetic structure in a mag- 
netic field along the c-axis, we measured the intensity of 
28 magnetic Bragg peaks in the (h, k, 0) plane for T=8K 
and H=5T. The data is best described by the basis vec- 
tors of the irrep T4, and the magnetic structure is given 
by 

rh£* = (1.60(4), 0.08(3), 0(2)) Ms 

m v c * = (0, 0.02(5), 0(2)) Ms (19) 



1. High-temperature incommensurate (HTI) structure 

For temperatures between Thl and Tpjj, Bragg reflec- 
tions occur at the (2n + 1 ± q, 2m + 1,0) and (2?i + 1 ± 
q, 2m + 1, 2m + 1) positions. The intensities of 170 mag- 
netic Bragg reflections were measured and can be ex- 
plained with a magnetic structure belonging to represen- 
tation T4 given in Table IVTTT1 The magnetic structure at 
T = 7K is given by 



rh^ 4 = 



(1.9(1), 0.2(1), 0.2(2)) Mfl 
(0, 0.0(1), 0.2(2)) jiB , 



(18) 




10 2 4 6 8 10 

T(K) 



The quality of the fit is given by i? p =0.14 and x 2 =6-7. 
The effect of the field along the c-axis is thus to reduce 
the ordered moment along the a-axis compared to the 
zero-field structure at a somewhat lower temperature. 
This may be because the field induces a uniform moment 
along the c-axis which reduces the amount of moment 
available along the a-axis. 




FIG. 11: As for Fig. M but for H = 5T. 



FIG. 12: Field dependence of magnetic Bragg peaks at 
T = 8.4K. Integrated intensities were obtained as for Fig. 
1101 Field dependence of the integrated intensity of (a) the 
AF (1,1,0) reflection, (b) the incommensurate (1 — q, 1,0) 
reflection and (d) the incommensurate (1 + q, 1, 0) reflection, 
(c) shows the field dependence of the incommensurate wave 
vector q. 
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FIG. 13: Temperature dependence of /i-integrated intensity 
of the AF (1, 1,0) reflection in a field of H = 8.5T, applied 
along the c-axis, showing the absence of a second order phase 
transition when cooling from the paramagnetic phase to the 
low-temperature commensurate phase. The data instead re- 
flect a cross over phenomenon for T ~ 8K. 




H(T) 



FIG. 14: Field dependence of effective staggered moment at 
three different temperatures, obtained from the (1,1,0) re- 
flection by taking the square root of the peak intensity. 



magnetic structure at 5K is given by 

= (1.6(1), 0.03(10), 0.01(7)) /z s 

rh^i = (0.0(5), 1.3(1), 0.1(1)) [lb 

m£« = (0, 1.4(1), -0.04(9)) ft B 

rh^i = (-2.2(1), 0,0) MB . (20) 

The quality of the fit is given by i? p =0.19 and x 2 =7. 
The corresponding spin structure is shown in Fig. 115b 
and consists of elliptical a — b plane spirals on spine and 
cross-tie sites. 

The structure at T = 5K thus consists of spirals on 
the spine and cross-tie sites, propagating along the a- 
axis with moments in the a-b-plane, as shown in Fig. 
I15H . An inspection of the spin structure suggests that 
nn interactions between Ni on adjacent spines are AF, 
both within and between Kagome planes. The moment 
in the c-direction is zero within the error bar, indicating 
the presence of a spin anisotropy which forces the spin 
into the ab plane. 

To determine the LTI magnetic structure in a mag- 
netic field along the c-axis, we measured the intensity of 
28 magnetic Bragg peaks in the (h, k, 0) plane for T=6K 
and H=2T. We obtained best agreement with the exper- 
imental data for basis vectors belonging to the irreps Ti 
and T^. The structure is given by 

n£« = (2.5(1), -0.1(2), 0(2)) /x B 

= (-0.5(3), 1.1(1), 0(2)) fi B 

n£« = (0,0.0(1), 0(2)) n B 

= (0.06(8), 0,0) mb, (21) 

The quality of the fit is given by i? p =0.29 and x 2 =9.3. 
Qualitatively, these parameters are similar to those of the 
zero-field structure. We remind the reader that Eq. i|ll|) 
implies a phase difference of 7r/2 between the a- and b- 
components on the spine sites. The high-field LTI phase 
thus consists of a spiral on the spine sites and no moment 
on the cross-ties, possibly because a transverse field has 
a strong effect on the cross-ties which are more weakly 
coupled antiferromagnetically. In contrast to the HTI 
phase, however, the ordered moment at non-zero field is 
higher moment than at zero field. 



2. Low-temperature incommensurate (LTI) structure 

For temperatures between Tlc and Thl, Bragg reflec- 
tions were observed at the (2n + 1 ± q, 2m +1,0) and 
(2n + 1 ± q, 2m + 1, 2m + 1) positions. These are the 
same Bragg peaks as observed in the HTI phase, but 
their relative intensity has changed, indicating that the 
spins undergo a spin reorientation at Thl- The present 
diffraction data are consistent with either a Ti + r 4 or a 
T 2 + T 4 structure. We choose the former spin structure, 
not because it has a smaller value of x 2 , but because it 
is consistent with the appearance of ferroelectricity, as is 
discussed in Sec. lIVCTl below. With that assumption the 



3. P and C phase 

The C phase has the same symmetry as the P phase, 
and for fields along the c-axis, there is no phase boundary 
between the high-field phase and the zero-field phase for 
Tec < T < Tlc- The symmetry of the C magnetic 
structure can thus be determined in a high magnetic field 
along the c-axis. We measured a set of magnetic Bragg 
peaks in the (h,k,0) plane at T = 0.1K and H = 8T, 
and we found that the data is best described by irrep T7 
and the parameters 

= (2.4(1), 0,0.0(5)) mb 
n£ T = (0,0.8(1), 0(1)) mb. (22) 
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The quality of the fit is given by P p =0.23 and \ 2 =17. 
When we allowed the fit to include also Ti parameters, 
the components of the magnetization were found to be 
0.1(l)/is and therefore statistically insignificant. 

A magnetic field along the c-axis induces magnetic or- 
der even in the paramagnetic phase. At T = 10K and 
H = 8T, the spin structure is best described by the irrep 
I\ with the following parameters: 



(0.61(2), 0, 0.0(2)) (xb 
(0, 0.37(3), 0.0(4)) /j,B 



(23) 



The quality of the fit is given by P p =0.31 and x 2 =2.49. 
Thus the field-induced magnetic structure at T = 10K is 
described by the same irrep as at T = 0.1K and H = 8T. 
This is corroborated by the absence of a phase transition 
upon cooling in a high field, as shown in the temperature 
dependence of the (1, 1, 0) Bragg peak shown in Fig. El 



Ni sublattice 




FIG. 15: Lattice (a) and magnetic structure (b-d) of various 
phases. The size of the + and — signs correspond to the com- 
ponents out of the page and into the page, respectively. In 
panel (d) the canting of the C phase is magnified for legibility. 
For the HTI and LTI phases, the a-component has a wave- 
length approximately 1.4a so that q ~ 0.27, as indicated in 
Fig. Q 



G. Ferroelectricity 

Recently NVO has been shown to have remarkable fer- 
roelectric behavior^ A spontaneous polarization P has 
been found to appear only in the LTI phase. In other 
words, upon cooling into the LTI phase, the ferroelec- 
tric order parameter develops simultaneously with the 
LTI magnetic order parameter. In Ref. it was pro- 
posed that the multi-component order parameter asso- 
ciated with this phase transition requires the trilincar 
coupling, V, where 



a 7 (J LTI (7 HTI + a-yOLTlOHX! 



P, 



(24) 



Here P 7 is the 7-component of the spontaneous polar- 
ization and the ct's are complex-valued order parameters 
which describe the incommensurate long range order as- 
sociated with irrep T4 for the HTI phase and with irrep 
Ti for the additional order parameter appearing in the 
LTI phase. Of course, V must be invariant under the 
symmetry operations of the paramagnetic phase^i*^ As 
we have seen, the magnetic order parameters satisfy Eq. 
(151). which here is 



I(7A = a\ 



(25) 



where A denotes either LTI or HTI. Using this relation, 
we see that the invariance of V under spatial inversion 
implies that et 7 is pure imaginary: a 7 = ir 7 , where r 7 is 
real valued. Thus we may write 



V 



2 ^2 r 7 P 7 |cr L Ticr H Ti| sin( 



'HTI 



<hn) , (26) 



where the phases are defined by a a = |o"A|e l ^ A . 

To be invariant under the operations of the little group 
P 7 must transform like I^ ® Ti. Referring to Table IVll 
we see that this means that P 7 has to be odd under 2 a 
and m ac . The first of these conditions means that r 7 
can only be nonzero for 7 = b or 7 = c. The second 
condition means that r 7 can only be nonzero for 7 = 6, 
in agreement with the observation^ that the spontaneous 
polarization only appears along the b-direction. Had we 
chosen the irrep T2 for the new LTI representation in 
Eq. H20|) . we would have incorrectly predicted the spon- 
taneous polarization to be along the c-direction. 



H. Magnetization and Susceptibility Measurements 

Although neutron diffraction enables one to fix many 
details of the magnetic structure, it is hard to obtain 
the bulk magnetization from these measurements because 
the signal associated with bulk magnetization is buried 
in the nuclear Bragg peaks. Accordingly, we summarize 
here the results for the zero wavevector magnetization, 
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FIG. 16: M versus H along the three crystallographic axes 
for a sequence of temperatures. 



M, measured with a SQUID magnetometer, as a func- 
tion of field for fields along each of the crystallographic 
directions shown in Figs. 1161 and HTI 

The following features are noteworthy. In Fig. one 
sees that for H along the c-axis, there is a range of tem- 
perature (corresponding to the C phases) in which M 
vs H does not extrapolate to zero for H — > 0. This is 
the best measurement of the weak FM moment in this 
phase. The C phase may also have a finite remnant 
magnetization though further measurements are needed 
there. From these data one sees confirmation of the 
phase boundaries obtained from specific heat measure- 
ments (shown in Fig. ^) which here are signaled by 
discontinuities in the magnetization as a function of H 
(when a phase boundary is crossed.) In Fig. El we show 
M/H, measured for a small field H = 0.1T. This quan- 
tity will be nearly equal to the zero field susceptibility 
except when it probes the spontaneous magnetization, 
as when H is along (0,0,1) and T < 4K and the system is 
in the C or C phase. It is remarkable that there are no 
visible anomalies associated with the phase transitions 
involving the HTI phase. 



V. THEORETICAL INTERPRETATIONS 

Here we discuss the simplest or "minimal" model that 
can explain the experimental results for NVO. In gen- 
eral, the higher the temperature the simpler the model 
needed to explain experimental results. Accordingly, we 
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FIG. 17: x f° r magnetic field applied along each of the crystal- 
lographic directions in Bohr magnetons per Ni ion per Tesla 
versus T, where \ is defined by \ = M{B. = O.YT)/H. 



discuss the theoretical models for the phases in the order 
of decreasing temperature. That way, we will see that the 
minimal model is constructed by sequentially including 
more terms as lower temperature phases are considered. 



A. The high temperature incommensurate (HTI) 
phase 

1. Competing nearest and next nearest neighbor 
interactions 

Incommensurate phases usually result from competing 
interactions. In the HTI phase we found spin ordering 
predominantly on the spine sites [see Eq. (|18|) ]. so the 
minimal model will deal only with the spine spins. Since 
the incommensurate wave vector is along the spine axis, it 
is reasonable to infer that this incommensurability arises 
from competition between nn and nnn interactions along 
a spine. We already indicated that such competition is 
plausible, in view of the nearly 90° angle of the nn Ni- 
O-Ni bonds. Thus the minimal model to describe the 
incommensurate phase is 



T^spinc — ~ ^ ^""^ Jna ^""^ S a (t) S a (jC + S n ) 



n—1 a 

7~tA + 7~L SS ' , 



r,<5„ 



(27) 



where J\ a and Jia represent the (a priori anisotropic) 
nn and nnn exchange interactions, a is summed over 
the axes a, b and c, r is summed over only spine sites, 
and S n = ±(rt/2)aa are the first and second neighbor 
displacement vectors for the spine sites (remember that 
there are two spine spins along each spine axis in the 
unit cell, so the nn distance is a/2). Also H ss ' represents 
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(probably weak) AF interactions between nn in adjacent 
spines: J oa for spins at a distance 6/2 in the b-direction 
and J ca for spins at a distance c/2 in the c-direction. 
(In this minimal model, for the purpose of Fourier trans- 
formation, we place the lattice sites on an orthorhombic 
lattice for which the a — c planes are not buckled, but 
we retain the same interactions between spins as for the 
buckled lattice. In this way our model gives the correct 
thermodynamics, even though it should not be used to 
obtain scattering intensities.) 

To have competition, we must have J 2a > (at least 
for the relevant a, which turns out to be a — a). Given 
that at low temperature we end up with a commensurate 
AF ordering of the a-spin components, it is also reason- 
able to assume that J\ a > (at least for a = a). If 
the spines did not interact with one another, they would 
form an array of independent one-dimensional systems 
for which thermal and/or quantum fluctuations would 
destroy long-range order. To understand the coupling 
between adjacent spines, note that a displacement (b/2) 
(from a spine in one Kagome layer to a spine in an ad- 
jacent Kagome layer) takes site #1 in one unit cell into 
site #4 in another unit cell and vice versa. Similarly this 
displacement takes site #2 in one unit cell into site #3 
in another unit cell and vice versa. From Table VII (or 
VIII) one see that for the active irreps #4 and #1 such a 
nearest neighbor displacement along b corresponds to a 
change in the signs of m sa and m so . The fact that these 
are the dominant order parameters for the spine sites 
therefore strongly suggests that the nn interspine inter- 
actions along b are antiferromagnetic. The evidence for 
antiferromagnetic interactions along c is almost as com- 
pelling. For these nn pairs (sites #1 and #2, or sites #3 
and #4) the largest order parameters of the LTI phase 
(see Eq. (PU)l )B. i. c. the x components of spine spins or- 
dered according to irrep #4 and the y component ordered 
according to irrep #1 are both antiferromagnetic along 
c. While it is true that the ordering of the x-component 
according to irrep #1 is ferromagnetic along c, this com- 
ponent of the order parameter is indistinguishable from 
zero. Thus we conclude that the nn interactions J\ a for 
all directions a are antiferromagnetic. 

Also, in Eq. (|27() Ha is a single ion anisotropy, 

Ha = -^T4 Q ]TS Q (r) 2 . (28) 

a r 

The continuum mean-field phase diagram of this Hamil- 
tonian is known, for both exchange and single ion 
anisotropics i^aS To determine which phase orders as the 
temperature is lowered from the paramagnetic phase, it 
is sufficient to look at the quadratic terms in the expan- 
sion of the free energy per spin in the Fourier components 
of the spins, 

F spi„e = 1 ^2J2[XsAP)r 1 S a (p)S a (-p) , (29) 

p a 

where 



is the a-component of the inverse susceptibility of the 
spine sites associated with the Fourier component 

S a (p) =Y,S a (r)e-^ r /M (31) 

r 

(with the sum over all the Af = 4iV uc spine spins in the 
lattice), while J a (p) is the Fourier transform of the ex- 
change interactions. In Eq. (|30|l . C is the Curie constant, 

C = S(S +l)/3 = 2/3. (32) 

(We measure energy in temperature units, which 
amounts to setting the Boltzmann constant k = 1.) For 
our simple model 

J Q (p) = 2[J la cos(p a a/2) + J 2a cos(p a a) 

+ J ba cos(p 6 6/2) + J ca cos(p c c/2)] . (33) 

As T is lowered, the first phase to order will involve the 
order parameter S a (p) for which [x s ,a(p)] _1 first van- 
ishes. For J\ a > 4| J2a\, this happens for 

p = Q = 27r(l/a,l/M/c), (34) 

implying a simple two sublattice antiferromagnet along 
each spine chain (with the orthorhombic unit cell con- 
taining four spins in each sublattice, and with the spins 
inside the unit cell varying with e l( 3°' r ). However, for 
J2a > I J\a |/4, the minimum in J Q (p) occurs at the in- 
commensurate wave vector 

p 0a = 2ir(q 0a /a, 1/6, 1/c). (35) 

The modulation wave vector qo a for the a spin compo- 
nent (in r.l.u.'s) is given by 

cos(irq 0a ) = -Ji a /(AJ 2a ) , (36) 

and thus the susceptibility for this wavevector becomes 

[Xs^Poa)]- 1 = T/C-2[J 2a + J? Q /(8Jaa) 

+ Jba +Jca + A a }. (37) 

Experimentally, we know that the leading order param- 
eter in the HTI phase concerns S a (poa)- In this phase, 
the neutron diffraction data also give a wave vector which 
varies slightly with temperature, close to qo a = 1 ± q ~ 
0.72 or 1.28. To avoid confusion, from now on we use 
the notation qo = 0.72. Using this approximate value, 
Eq. (123) gives 

J la » 2.55J 2a . (38) 
Thus we conclude that at this mean field level one has 

T PH = 2C[J 2a + Jl 2 a /(8J 2a ) + Jba + Jca + A a ], (39) 

and we end up with a longitudinally modulated spin 
structure, with S'h(r) = S c (r) = and, along a single 
spine, 



[X^a(p)]" 1 =T/C + J a (p)-2A, 



(30) 



S a (x) = Scos(2Trq x/a + <j) n ), (40) 
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where the phase <f> n may depend on the spine index n. 
The n-dependence of <f> n would be determined by the 
inter-spine coupling. For our model the neutron diffrac- 
tion data indicates that adjacent spines are antiferromag- 
netically arranged, so that for a general spine site we may 
write 

S„(r) = S , cos(po-r + o ) , (41) 

where the transverse components of po are fixed as in Eq. 
(|33|l and we adopt the notation p = po a - When one goes 
beyond continuum mean-field theory, it is found that in- 
stead of go being a continuous function of J la and J 2o one 
obtains a 'devil's staircase' dependence of the wave vec- 
tor on the control parameters ^SiSi This treatment also 
shows that <fio can not be fixed arbitrarily (as it would be 
in continuum mean-field theory). Furthermore, critical 
fluctuations will reduce the actual Tph by a factor which 
depends on the spatial anisotropy. 

2. Dzyaloshinskii-Moriya interactions 

In addition to the leading a-component of the spine 
spins, Eq. (|TB|) also indicates a small, but non-negligible, 
spine moments along the b-axis and along the c-axis. 
Such moments follow directly from the Dzyaloshinskii- 
Moriya (DM) interactions^^ between spine spins. For 
simplicity, we consider here only the nn DM interactions 
along the spine, 

W DM = ^D(r).S(r)xS(r + fl a/2), (42) 

r 

and the symmetry of the lattice dictates that the DM 
vectors behave as 

D(r) = (0,Ae iQ °- r ,^e iPe - r ), (43) 

where Qo is the wave vector for the two-sublattice com- 
mensurate wave vector [Eq. (J331], while P c = (0, 0, 2tt/c) 
represents AF ordering along the c-axis. Next nearest 
neighbor DM interactions, discussed in II, do not change 
the qualitative results presented here. The Fourier rep- 
resentation of TLdm yields a free energy of the form 

F DM = J2[ D bCos( Pa a/2)S a (p)S c (-p - Q ) 
p 

-iD c $m(p a a/2)[S a (p)S b (-p-P c )]+h.c. (44) 

Introducing quadratic terms in the transverse spine 
spin components, as in Eq. 1|29[1 . we can now minimize 
the free energy and find these transverse components: a 
non-zero 5* a (po) generates 

<S c (Po + Qo) « -2D b Xs,c(po + Qo) cos(7rqo)S 

S&(Po + Pc) « 2iD cXs APa + Pc) sin(7rg )5 a (po){45) 

consistent with the signs and phases of fh b s in Table IVTTT1 
In fact, for all the group representations, all the internal 



structure within the orthorhombic unit cell can be repro- 
duced using factors like e lpr , with the three wave vectors 
Po, Po + Qo and po + P c (with possibly different po's for 
different representations). 

Using the values from Eq. I|18|) , and p a = 2nq§ / a with 
q « 0.72, Eq. yields 

D c « -0.07[ X ^(po+Pc)r\ 

D b « 0.04[ Xs;C (p + Qo)]" 1 . (46) 

Near T PH , [x s , a (Po)] _1 is small, but [Xs,b(Po + Pc)] -1 
and [Xs,c(Po + Qo)] 1 remain finite. Since Tph = 9. IK, 
it is reasonable to assume that all these inverse suscep- 
tibilities are of order 10K. Thus we have the estimates, 
D b ~ 0.4K and D c ~ — 0.7K. Below we conjecture a set of 
exchange parameters, which give [xb(Po + Pc)] 1 ~ 13K 
and [xc(Po + Qo)] -1 ~ 47K, changing these estimates 
into D b ~ 0.5K and D c 3K. 

3. The spine-cross-tie interactions 

The representation also allows some small incom- 
mensurate moments along the b and c axes on the cross- 
tie sites. Indeed, the experimental values in Eq. (|18fl also 
allow for such moments, alas with large error bars. Here 
we discuss the possible theoretical origin for these mo- 
ments, namely the anisotropic spin interactions between 
the spine and the cross-tie spins. In all the phases, the 
two spine chains which are nearest neighbors to a given 
row of cross-tie sites in an a — c plane (in this section 
we ignore the buckling, which does not affect the present 
considerations) have antiparallel a-components of spins. 
Thus even in the incommensurate phases, the isotropic 
nn spine-cross-tie interaction is frustrated, and one needs 
to add (symmetric and antisymmetric) anisotropic spine- 
cross-tie interactions. In the simplest approach, these 
interactions can be written in terms of an effective in- 
ternal field produced by the four spine spins (si to s4) 
surrounding a cross-tie spin c2 (see Fig. |2J, so that 

H sc = -J2KS c a (c2), (47) 

a 

with 

s4 
P si 

where the matrix J(i) (which relates to the coupling be- 
tween the spins at si and at c2) contains both symmet- 
ric (pseudo-dipolar, PD) and antisymmetric (DM) off- 
diagonal terms, whose signs depend on i (see II). Ignor- 
ing the interactions among cross-tie spins, each such spin 
will follow its local field, 

^(c2) = Xc, Q ^, (49) 

where Xc,a is the a-component of the cross-tie suscepti- 
bility. 
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The above analysis yields the spin components on each 
cross-tie site in terms of the four surrounding spine spins. 
Taking also into account the variation of the matrices 
J{i) for different plaquettes, and assuming only linear 
response for the cross-tie spins, we end up with 

^(Po) = -4x C)C (j ac - 4)5 a (p )sin(7rq /2), 
5* fe c (p + P c ) = 4x c ,6(iafc + 4)5* a (po)sin(7r 9 o/2) ,(50) 

where j a p and c? 7 are the symmetric (PD) and anti- 
symmetric (DM) elements of the matrix J{\). Using 
the experimental values of the cross-tie moments from 
Eq. (|18(l . we end up with 

(jab + d c )Xc,b = -0.06 ±0.06, 

Uac - d b ) Xc ,c = ± 0.06 . (51) 

Ignoring the very weak interaction among the cross-tie 
spins, we can use the free spin Curie susceptibility, Xc,a ~ 
1/T ~ 0.2/K, hence {j ab + d c ) k -0.3(3)K, ( Jac ~ d b ) sa 
0.0(3). However, the uncertainties of these estimates arc 
very large. 

B. The low temperature incommensurate (LTI) 
phase 

Mean- field theory 2 ^^ also indicates that, if the uniax- 
ial anisotropy energy is not too large, then as the temper- 
ature is reduced further below Tpjj, a second continuous 
phase transition occurs at Tjjl, below which transverse 
modulated order appears, leading to an incommensu- 
rate state with elliptical polarization. Again, this seems 
to be consistent with the experiments: Eq. I|20|l indi- 
cates a growing ^-component on the spine spins and a 
growing a-component on the cross-tie spins. For sim- 
plicity, we again concentrate only on the spine spins. 
One can understand this transition intuitively by real- 
izing that as the temperature is lowered the fixed length 
constraint on spins (embodied in the quartic terms in 
the Landau expansion) becomes progressively more im- 
portant. Since the experiment indicates transverse or- 
dering in the b-direction [see Eq. H20|) ]. we conclude 
that the anisotropies still require that S c (r) = 0, and 
we write a Landau expansion in both S a (p) and Sf,(p). 
A priori, the quadratic terms in Sa(p)| 2 and in |Sb(p)| 2 
could have minima at slightly different wave vectors pi 
and p 2 respectively [see Eq. (J3BJ)]. The quartic cou- 
pling between the two spin components would then be 
of the form |S' a (pi)| 2 |S'b(p2)| 2 - Depending on the ratio 
of the amplitude of this term to those of S^pi)! 4 and 
5 f b(P2)| 4 , one could end up with a second transition, at 
Tjjl < Tph, at which Sb (P2) would begin to order (this 
is similar to the usual appearance of a tetracritical point, 
see e.g. Ref . I30I) . However, the quartic terms also contain 
a term of the form [S a (pi) 2 Sb(— P2) 2 + h.c. ]<5(pi — P2) = 
2|S , a (pi)| 2 |5 , b (p 2 )| 2 cos[2(0 Q - 4>b)]8{vi - p 2 ), where we 
have written S a (p) = \S a (p)\e t ^ a . Usually, the coeffi- 
cient in front of this term is positive, and therefore it is 



minimized when cos[2(</> a — 4>b)] — —1, and it lowers the 
free energy only if pi = p 2 . If the values of pi and P2 
are sufficiently close to one another, then the presence 
of this term results in locking the wave vectors to each 
other: Pi = P2 = Po ; consistent with the experiment. 

Thus, within this theory we expect spin ordering of the 
type 

S a (r) = S a cos(p ■ r + <j) a ), 

Sf,(r) = 5 h cos(p • r + fa) , (52) 

with S b growing continuously from zero as the tempera- 
ture is lowered through Thl, and with 2((f> a — 4>b) an odd 
integer multiple of 7r. The fact that the two order pa- 
rameters are out of phase can be understood intuitively: 
if longitudinal and transverse spin components are com- 
bined, they are closer to obeying the fixed length con- 
straint if they are out of phase with one another. 

Within the Landau theory, we would expect \S a \ 2 to 
grow at (Tph — T) down to X"hl- Below that temperature, 
1 5b 1 2 grows linearly with (Thl — T), while the slope of 
IS'al 2 versus T decreases^ Qualitatively, this is what one 
observes in Figs. 151 1101 and II II 

The above theory ignores higher harmonics of the in- 
commensurate order parameters. Although these har- 
monics may be negligible close to Tph, they may grow 
at lower temperatures. Therefore, we went beyond the 
Landau expansion and numerically minimized the mean 
field moments on long chains of spin- 1 ions with the inter- 
action 12711 . with an isotropic exchange which has J\ — 
2.55J2 and with a uniaxial single-ion energy A a = K 
(see II for details). The resulting phase diagram, shown 
in Fig. 1181 is qualitatively consistent with the above pic- 
ture. It is intuitively clear that the range in tempera- 
ture over which the LTI phase is stable decreases as the 
anisotropy is increased. (Large anisotropy disfavors the 
existence of transverse spin components.) 

Having derived the leading order parameter, we can 
now evaluate secondary spin components, arising due to 
either the spine-spine DM interaction or the spine-cross- 
tie PD interactions. Basically, such an analysis will gen- 
erate all the other spin components that are allowed by 
representations T4 and Y\. In the above discussion, we 
represented T4 by its largest component 5 Q (po) = Tn%£, 
and Ti by Sh(po) = m s b- Once this is done, one could 
in principle deduce more information on the PD and DM 
coupling constants. Unfortunately, the only data avail- 
able in the LTI phase is at 5K, and these data show ap- 
parently large values of the cross-tie spins. These large 
values are surely beyond our linear response treatment, 
and therefore we are not able to use them for identify- 
ing the coupling constants. At the moment, we do not 
understand this apparent fast saturation of the cross-tie 
spins. 
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FIG. 18: Mean-field phase diagram for the isotropic J\-Ji 
model with easy-axis anisotropy scaled by K = A a - Here 
"HTI" denotes a longitudinally polarized incommensurate 
phase and "LTI" an elliptically polarized incommensurate 
phase. In both phases the modulation vector is given by Eq. 

"AF" denotes a two-sublattice collinear AF phase. For 
large anisotropy this model reduces to the ANNNI modeli— 



C. The commensurate antiferromagnetic C phase 



Structure and symmetry 



As noted in Refs. and as confirmed in our mean 

field analysis (Fig. I18fl . lowering the temperature yields 
a first order transition between the LTI phase and a 
commensurate phase. In this commensurate phase, the 
spine spins return to be along the a-axis (like in the HTI 
phase), but now they all have the same length. At zero 
field we identify this phase with the C phase. At finite 
field along the c, this phase coincides with the P phase. 
For our minimal model 7i S pinc, the order parameter in the 
P phase is the simple two sublattice staggered moment, 
N s , a = S a (Q ). 

In addition to gaining energy from the fixed length 
constraint, one also gains energy from the appearance 
of a weak FM moment. This moment arises due to the 
anisotropic PD and/or DM interactions. We have already 
noted that the symmetry of the P phase is confirmed by 
the lack of a phase boundary to the paramagnetic phase 
at high magnetic field along the c direction. (In a mag- 
netic field along c, the paramagnetic phase must be in- 
variant under a two-fold rotation about the c-axis and 
under inversion because H is an axial vector. Thus the 
magnetic structure must be invariant under these oper- 
ations but should change sign under two-fold rotations 
about the a or b-axis. This means that in the paramag- 



netic phase with a field applied along the c axis we must 
have only the representation Tj of Table IVl) The above 
conclusions are supported by the structure observed at 
T = 0.1K and H = 8T. 

A similar symmetry analysis shows that at high tem- 
perature we have only the irrep of Table Ivl for H along 
a and for H along b. Thus, for H along a or b, a 
comensurate phase containing Tj can only be reached by 
crossing a phase boundary. This is also consistent with 
the data: for fields along the a- or b-axis, the paramag- 
netic phase and the C phase are always separated cither 
by the HTI or by both the HTI and LTI phases. 

As stated, an external uniform field H along c gen- 
erates a nonzero magnetization on both the spine and 
cross-tie sites. The DM and PD interactions then gen- 
erate a nonzero staggered moment on the spine sites, 
N St a. In the paramagnetic phase, this moment will be 
approximately linear in the field. This explains the 
nonzero intensity at Q = (110) in the P phase, seen in 
Figs. ^2 EH! an d El Within the Landau theory, sym- 
metry allows a biquadratic coupling of this staggered 
moment to the incommensurate order parameters, e.g. 
iV? i0 |5 (po)| 2 - Near the P-HTI transition, this term 
renormalizes [Xa(Po)] -1 by an amount of order H 2 , thus 
shifting Tpu by such an amount. 

At higher H , the coupling of the incommensurate HTI 
order parameter to both M and N also renormalizes the 
quartic term ^(po)! 4 , yielding a tricritical point when 
this term turns negative. 

Below T PH , |5"a(po)| 2 grows as (T PH - T). The above 
quartic term then induces a corresponding change in the 
inverse staggered susceptibility [Xa(Qo)]^ 1 , causing the 
change in slope of 7V Sj0 versus T seen in Fig. EJ More 
of these calculations are given in II. 



2. Weak ferromagnetic moment 

Now we discuss the possible sources of the small 
FM moment observed in the C phase and possi- 
bly also the C phase. This can arise either from 
the DM interactionsS&Sl or from the pseudo-dipolar 
interactioniiSiSi^SiS We start by considering the DM in- 
teraction. Equation (|45|l (with go replaced by 1) shows 
that when po is replaced by Qo then the spine staggered 
moment N SM generates a FM moment 



2D bXs ,c{°) N s,*- 



(53) 



(From now on we distinguish between spine and cross-tie 
properties by the subscripts s and c). 

Using our rough estimates D\, ~ 0.4K and Xs.c ~ 0.1K, 
and ignoring the spine-cross-tie coupling (see below), this 



gives M s 



0.08iV, 



0.12/j.b- This spine value is 



within the error bar given for this moment in Eq. 12211 . 
The same equation also allows for a similar moment on 
the cross-tie sites. That same equation then also implies 



a FM on the cross-ties of order M r 



-0.1/ifl. 
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Similarly, the off-diagonal exchange spine-cross-tie in- 
teraction will generate a FM moment in the c-direction 
on the cross-tie sites: Eq. (l50"|) implies 



Mr 



-4x c ,c(ioc - d b )N s 



(54) 



and the estimates (j ac — <4) = 0± 0.3K and Xs,c ~ O.IK 
give M c c = ± 0.2[1b, which are consistent with the 
above estimate. 

These estimates are justified only if we ignore the di- 
rect exchange coupling between the FM moments M s and 
M c . Although this coupling cancelled out for the AF and 
incommensurate moments, the diagonal elements in the 
matrices J{i) do add up for the FM moments. Assum- 
ing isotropic diagonal elements, this generates a coupling 
4 J av M s -M c per cross-tie site, where J av is the spine-cross 
tie nn interaction. A complete analysis then requires a 
minimization of the free energy with respect to the two 
FM moments. The situation here is very similar to that 
which occurred in S^C^C^O?^ which also has a weak 
FM moment in an AFM phase and involves two types of 
magnetic ions (copper). Following Ref. |33, we assume 
that at low field in the C phase the staggered moment on 
the spine sites is practically saturated. The free energy 
of the small FM moments per Ni ion can then be written 
as 



3F 



Mr. 



1 Xs 



2Xc 



2H S • M s - H c • Mc, (55) 



where we denote J pG \ = j ac — db and 
U s = U-2D b N s , a c, 

H c = H - 4J av M s - 4J pc \N Si aC. (56) 

Minimization with respect to the magnetizations then 
yields the total magnetization per Ni ion, 



M a = M S a , c + XaHa, 

with the average zero-field moment 

x(4J av JpdX c , c - Db) - JpdXc 
and susceptibility 

^Xs,a ~t~ Xc,a 8<^avXs,aXc,a 



(57) 



3 (! - &JivXs,aXc,a) 



(58) 



(59) 



Note that For J av = 0, this reduces to the trivial M a = 
(2/3)M s , Q + (l/3)M c , a . 

Experimentally, we seem to observe a weak tempera- 
ture dependence of Mq, with a stronger T-dependence in 
Xs,a- In the C phase, the spine spins are ordered anti- 
ferromagnetically, and therefore the Xs.a's have a weak 
T-dependence (both for the transverse and the longitudi- 
nal susceptibilities). Therefore, the temperature depen- 
dence in the above expressions probably comes mainly 



from Xc,a- Expanding Mq to leading order in Xc,c, we 
note that this order vanishes if J p d = 2x s ,cJavDt- In the 
AF C phase, Xs,c is the transverse susceptibility, which 
is of order l/[4(Ji - J 2 ) + 2(J b + J c + A a )]. It is also 
reasonable to assume that J av is of order Ji, given the 
similar geometry of the nn ss and sc bonds. This would 
imply that J av Xs.c is of order one. Since we also found 
that D b and J p d are of the same order of magnitude, the 
above relation may in fact hold approximately, and this 
could explain the weak T-dependence of Mq. However, 
this is all highly speculative at this stage. 



D. Phase boundaries in the H-T plane 

1. Qualitative discussion 

Now that we have identified the magnetic structure 
of all the ordered phases, we can give a qualitative dis- 
cussion of the various boundaries in the phase diagram. 
In the various incommensurate phases, the free energy 
must be an even function of the magnetic field. Since 
the field generates a uniform magnetic moment, it com- 
petes with the incommensurate order. Therefore, we ex- 
pect the lines between the high temperature paramag- 
netic phase and the HTI phase to behave at low field like 
Tpn(H) k Tq—AH 2 , with the coefficient A depending on 
details. As stated above, this coefficient arises due to the 
biquadratic coupling of the HTI order parameter to N 2 a 
and to M 2 . The data seem consistent with this, showing 
a small value of A. The observed phase boundaries be- 
tween HTI and LTI also seem almost field independent, 
remaining roughly at a constant distance below Tpn(H). 
To leading order, the instability which yields the trans- 
verse incommensurate order is mostly dominated by the 
size of the longitudinal order parameter, which depends 
mainly on this temperature difference. 

We next discuss the (first order) boundaries between 
the phase C and the two incommensurate phases. We 
start with the two upper panels in Fig. ^ where the field 
is not along the weak FM moment in phase C. A guid- 
ing principle in this discussion is that the Zeeman energy 
of an antiferromagnet (or a general incommensurate or- 
dered state) in a uniform field, — (l/2)xH 2 , is larger in 
magnitude when the field is perpendicular to the stag- 
gered moment than when it is parallel to that moment 
(because x± > X\\)- Now consider the H-dependence of 
the transition temperature Tlc from the LTI phase to 
the C phase. When the field is along (100), the com- 
mensurate phase gains very little Zeeman energy (due to 
X\\), whereas for the incommensurate phase, the field is 
partly transverse, because the spins are in a cone. Indeed, 
Fig.QJi shows a transition from C to LTI as H increases, 
with a phase boundary which might be parabolic. In 
contrast, when the field is along (010), the commensu- 
rate phase gains the entire energy — (l/2)x±H 2 , whereas 
in the incommensurate phase, to the extent that the spins 
are partially in the b-direction, the lowering of energy is 
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less than that in the commensurate phase. We thus ex- 
pect a transition from the LTI phase to the C phase as H 
along (010) is increased, consistent with Fig.^i. Again, 
the phase boundary seems parabolic. 

In contrast, when the field is along (001), the com- 
mensurate phase is strongly favored because in addition 
to the energy involving the transverse susceptibility, the 
external field also couples linearly to the spontaneous 
FM moment, which is along the c-direction. This ex- 
plains why the C-LTI and the C-HTI phase boundaries 
are practically linear in the H-T plane only for H along 
(001). 

As discussed above, the smearing of the transition be- 
tween the paramagnetic phase and the C phase results 
from a bilinear coupling between the magnetization M c 
and the AF order parameter N s , a , implying that a uni- 
form field along (001) also acts as a staggered field on 
N Sja . Had it not been for this coupling, then the com- 
petition between HTI and C would result in a bicritical 
phase diagram, with the two second order lines (P-HTI 
and P-C) meeting the first order HTI-C line at a bicrit- 
ical point. Indeed, Fig.^J; does show the typical inward 
curvatures of the former two lines, typical for such phase 
diagrams^ 

2. Field along (001 ) 

We now make these arguments more quantitative. We 
start with the LTI-C and HTI-C phase boundaries. Since 
these represent first order transitions, it suffices to com- 
pare the free energies in the respective phases. Since the 
susceptibilities may have different values in the different 
phases, we use the superscripts C and / for the commen- 
surate and the incommensurate phases. We start with 
the field along (001). 

In the commensurate C phase, with a field along (001), 
we write the energy per Ni ion as 

Ec=E° c - \x% c {^H - H s f - \xcAH - H c ) 2 , 

= E° C - \x C c H 2 - HM - \xZ c H 2 s - \xccH 2 c ,(m) 

where = l-2J av x c , a , x£ Q = x£«/(l - 8 J% v X?, a Xc, a ) , 
H s = (2D b - 8J av J pd Xc,c)iV s , a) He = -4J pd W s ,a, M = 
(2/j, c x^ c H s + Xc, c H c )/3. Also, the total susceptibility 

is X? — ^(2m 2 X^ c + Xc,c)- It is easy to check that 
these expressions agree with those in Sec. IV C 21 At 
low temperatures, our isotropic Ji — J2 — A model yields 
Eq — |(— Ji + J2 — A), where A = A a is the coefficient 

of the single ion anisotropy, Eq. |2HJ>, and x% c — X?j_ ^ s 
temperature independent. 

In the incommensurate phase LTI, we may have dif- 
ferent values for Xs.a, which we now denote xl a- As 
we don't expect Xc.a to change, we now have the total 
susceptibility as 

xl = (ZfJ-ixl,* + Xc,a)/3 



_ ^X s .q + Xc.a &Ja,vXs,uXc,a . . 

3(1 - 8Jlx{ a Xc, a ) ' 1 j 

Since in this phase, N s<a is replaced by an oscillating 
term, the "spontaneous" moment on the spine spins, Mo, 
will be replaced by an oscillating term which will give no 
average contribution to the energy. For the same reason, 
we need to replace H 2 and H 2 by their averages over 
these oscillations, which we denote by jH 2 and SH 2 . 
Thus we write 

Ei = E a j - \x\H 2 - ^xl c H 2 - \5xccH 2 c . (62) 

For our J\ — J2 — A model, at low temperatures, we have 
E I = l[-J 2 - J!/(8J$)-A/2}. 

The first order transition between these two phases 
will occur when Eq = £7- Experimentally, it turns out 
that the difference Xc ~ xi is very small. Neglecting this 
difference, remembering that the experimental Mo seems 
practically temperature independent, and neglecting any 
temperature dependence of x't c an( i xl c > we obtain the 
transition field Hc.i at temperature T as 

M H C j = (x^c(Ti) - x c U T ) 

- 7Rlo(Ti)-X s yT)])i/ 8 2 /3 

+ (l-5)i? c 2 (xc,c(Ti)-Xc,c(T))/6, (63) 

where T\ « 3.8K is the transition temperature at H = 0. 
The apparent linearity of this transition line thus fol- 
lows from an approximate linear decrease of Xc,c(T) in 
the range 3.8K< T <6K, which also results in a linear 
dependence of x£ c (T) and xl,c( T )- 

To estimate the slope of Hci(T), we need detailed in- 
formation on the temperature dependence of the various 
individual susceptibilities. Experimentally, we have only 
information on the total susceptibility, Xc ~ xl ~ Xc ■ 
Therefore, we are not able to make any quantitative pre- 
dictions based on the observed slope. 

3. Other C-LTI phase boundaries 

For fields along (010) or (100), there is no linear term 
HMq in the analog of Eq. I|60l) . We also need to replace 
the relevant susceptibilities, according to the direction 
of the field. For a field in the a-direction, the phase 
boundary is thus given by 

\{x C a -xiWci 
= (x^m) - X C sAT) - lixljTx) xUT)])h 2 /3 
+ {l-5)H 2 (x c ATi)-XcAT))/a>, (64) 

yielding a parabolic dependence of T on the transition 
field Hci- As seen from Fig. El (Xa ~ Xa) IS larger and 
of opposite sign to (xf — x{)i explaining the shapes of 
the two other phase boundaries. 
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E. Temperature and orientation dependence of the 
susceptibility 

We now return to Fig. El To a good approximation, 
all the curves in this figure can be described by straight 
lines, as predicted in Eq. i|57|). The intercept M is zero, 
except for the high field data with H || (001). In the 
latter case, Mo extrapolates to a value which seems to be 
temperature independent. 

The data in Fig. 1161 exhibit two types of transitions: 
for the field along (100) (upper panel), there is a tran- 
sition from the C phase (low H) to the LTI phase (high 
H ) . The only change at this transition is a discontinuous 
increase in Xioo- This agrees with our expectations: as- 
suming that the change comes mainly from the ordered 
spine spins, this susceptibility is longitudinal in the C 
phase, and has some contributions from the transverse 
components in the LTI phase. In both phases, the spon- 
taneous moment Mo along (100) vanishes. 

For a field along (001) (lowest panel), there is a transi- 
tion from the LTI phase (low H) to the C phase, where, 
Xooi remains almost -ff-independent, but Mq exhibits a 
jump. The main gain in energy comes from the spon- 
taneous uniform magnetization. In addition, in the C 
phase, Xooi is fully transverse, whereas in the LTI phase 
it contains at least some longitudinal component. How- 
ever, this difference seems too small to be observed. 

The susceptibilities Xa all depend on temperature, 
with xoio and xooi decreasing as T increases, while xioo 
increases with increasing T. To analyze these results 
quantitatively, we must use Eq. (|59(l . which requires 
some assumptions on the separate susceptibilities on the 
spine and on the cross-tie sites. Qualitatively, it is rea- 
sonable to guess that the increase in x with decreasing T 
for H along (001) and (010) comes from the increase in 
the corresponding (possibly Curie-like) susceptibilities of 
the cross-tie spins, which are disordered. For fields along 
(100), the observed decrease in x with decreasing T is 
probably due to the strong decrease in the longitudinal 
spine susceptibility. However, a full quantitative analysis 
requires more experimental information than is currently 
available. 



F. Conjectured parameters 

The most robust conclusion from our experiments and 
analysis is Eq. (|38|) . for the ratio Ji a /J 2 a- Additional 
information concerning these exchange energies may be 
obtained from the high temperature susceptibility ten- 
sor, whose components are given asymptotically at high 
temperature as 

X aa ~C a /(T + O a ) , (65) 

where a labels the Cartesian component, (xa/3 must van- 
ish for a f3 for Cmca crystal structure.) For NVO the 



Hamiltonian is of the form 

h = - y^y^gq(i)/iB»s , a (?)ij r a 

+ \ M a7 (k,l)S a (k)S 7 (l) , (66) 

where M Q7 (i,z) is symmetric. Generalizing Eq. (|32[1 . we 
findi£ that 

C a = [(mb) 2 /3][8. 9q (s) 2 +4 5q (c) 2 ] , (67) 
where s denotes spine and c cross-tie, and 

4/Lt 2 6 

C a Q a = —^-^2^2g a (i)g a (j)M aa (ij) 

8=1 j^i 

„2 

+ -^-^2g a (i) 2 m aa (i 1 i) , (68) 

i 

where the sums count all the interactions in one primitive 
unit cell. If g(s) = g(c) (g « 2.3 34 ), then 

C a = ^ 2 B gl (69) 

and 

1 6 1 

6 « = gSS M «x(* s j')+ jgJ^Mcafti) • (70) 

i=l jjti i 

If, for instance, the exchange interactions are isotropic, 
then 

8 16 2 

S„ = -(Ji + J 2 + J b + Jo) + y^av - , (71) 

and the other components are 

Qb.c = ^(Jl + J 2 + Jb + Jc) + Y Jav + , (72) 

where J av is the isotropic nn spine-cross-tie interactions 
and the anisotropy constant K is defined by 

H A (i) = -K(Sf x - 2/3) 

K 2 2 2 

= ~ irPSix — &iy ~ Siz)- (73) 

Here we only attributed anisotropy to the spines. Now we 
consider what values can be obtained from the data. For 
high temperatures, a fit of the susceptibilities to Eq. I|65|l , 
as shown in Fig. HI yields @ aAc = 17±2K, 20±2K, and 
19 ± 2K, respectively. These values are not really consis- 
tent with an easy axis anisotropy. Rather they indicate 
an easy plane anisotropy with a smaller anisotropy which 
favors a over c. Using just the average value of 19K, we 
obtain the constraint 

8 16 

-(Ji + J 2 + J b + J c ) + y Jav wl9K. (74) 
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FIG. 19: Determination of the Curie- Weiss temperatures Q 
from the high-temperature susceptibility datai— 



We next estimate the anisotropy K. Figure ITH1 implies 
that the LTI phase appears between the HTI and the C 
phases only when K/Ji ~ 1/2. This value is also roughly 
consistent with the anisotropy of the susceptibility. 

Given these constraints, we propose a possible set of 
exchange parameters: 



Ji = 10, J 2 = 4, J av = 3 
Jb = J c = 2 , and K = 5 , 



(75) 



all in units of K. The uncertainties are probably of or- 
der 50%. Substituting these values into Eq. (|39|) then 
yields the mean field transition temperature Tpn ~ 20K, 
which is about 2.2 times the actual transition tempera- 
ture. This reduction from the mean field estimate must 
result from fluctuations. In ideal Kagome geometry one 
might think that J av and J\ would be identical. How- 
ever, in NVO the Ni-O-Ni bond angles^ for J x (95.0° 
and 90.4°) are further from 90° than those for J av (90.3° 
and 91.5°), which explains why J\ is larger than J av . Our 
estimates for the DM and PD parameters were already 
given in Sees. IV A 21 and IV A 31 



VI. CONCLUSIONS 



addition, for the LTI phase neutron diffraction data, on 
its own, does not unambiguously identify whether the ad- 
ditional representation which is characteristic of the LTI 
phase is Y\ or Y-i- However, the symmetry analysis of 
the spontaneous polarization (see Ref. |8() indicates that 
the correct choice is T\. 

2) We showed that a model having nearly isotropic in- 
teractions between nn and nnn on a spine and with sin- 
gle ion anisotropy can explain qualitatively the observed 
structures. (See Fig. 1181 ) In particular, the experimen- 
tal determination of the incommensurate wavevector ac- 
curately determines the ratio J1/J2 of nearest- to next- 
nearest-neighbor interactions on a spine. 

3) Using the experimental values for the P to HTI tran- 
sition temperature and the Curie- Weiss temperatures (of 
the high temperature uniform susceptibility), we obtain 
the estimates J\ ~ 10K and J 2 ~ 4K. 

4) From the appearance of small off-axis components of 
the magnetization in the HTI phase we obtain estimates 
in the range \D\ ~ 0.5K for the Dzyaloshinskii-Moriya 
interactions between adjacent Ni spine spins. 

5) From the appearance of a weak FM moment in the 
commensurate AF phases (and the appearance of trans- 
verse spin components in the HTI phase) we conclude the 
existence of either anisotropic symmetric exchange in- 
teractions or, more probably, Dzyaloshinskii-Moriya an- 
tisymmetric spine-cross-tie exchange interactions. Both 
interactions are permitted by crystal symmetry and we 
give very crude estimates of their values. 

6) We have provided a qualitative explanation for the 
shape of the phase boundaries between the various or- 
dered and paramagnetic phases. In principle the shape of 
these phase boundaries can be used to deduce additional 
microscopic interactions. However, there are currently 
too many unknown parameters to allow this program to 
be carried out. 

This work suggests several fruitful lines of future re- 
search, some of which are ongoing. Since the nn and 
nnn interactions along the spine are of the same order of 
magnitude, it would be interesting to alter the geometry 
of the coordinating oxygen ions. This might be done by 
applying pressure, or possibly uniaxial stress. Changing 
the ratio of these interactions would profoundly affect 
the magnetic properties of NVO. It is even possible that 
the magneto-electric properties could be grossly affected. 
Such experiments could greatly expand the emerging un- 
derstanding of the magnetoelectric behavior of frustrated 
quantum magnets. 



In this paper we have presented a comprehensive inves- 
tigation of the magnetic phase diagram of NVO. Here we 
summarize our conclusions and call attention to a num- 
ber of topics for future research. Our conclusions are 

1) The magnetic structures of the HTI, LTI and C 
magnetic phases of NVO have been closely determined. 
However, there are still some uncertainties in the struc- 
ture: some of the complex phases of the various complex 
order parameters are not determined with precision. In 
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TABLE IV: Irreducible representation of the group G v for the 
commensurate magnetic structure with v = (0,0,0). 



APPENDIX A: MAGNETIC GROUP THEORY 
ANALYSIS 

Using the space group symmetry of the crystal struc- 
ture we identify allowed basis vectors for a magnetic 
structure with the observed wave-vectors v. This is done 
by determining the irreps and eigenvectors of the little 
group G v of symmetry operations that leave the wave- 
vectors v invariant. The eigenvectors </> A of the A-th ir- 
reps r A were determined using the projector method^ 
They are given by 

A =£x A (s)ffW>), (Al) 
g 

where g is an element of the little group and (f> is any 
vector of the order parameter space. x A (f) is character 
of symmetry element g in the representation T A . 

The symmetry elements of the Cmca space group of 
NVO are given in Table [H] Note that this space group 
is nonsymmorphic because some of the group elements 
{C|a} consist of a reflection or rotation O and require a 
translation a equal to half a direct lattice vector. 



1. Commensurate structure 

The ordering wave-vector v = (0, 0, 0) is invariant un- 
der all these operations so that the little group G v con- 
tains all the elements of the space group. The group 
consists of 8 different classes and therefore has 8 irreps, 
all of which are one dimensional. We determined the 
classes and the character table of this group and these 
are shown in Table Hvl 

The eigenvectors were calculated using the projector 
method (See Eq. HAlfl .) and are given in Table Ivl 



2. Incommensurate structure 

For the ordering wave- vector v = (q,0, 0), the little 
group G v contains the following four elements of the 
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TABLE V: Irreducible representations for the commensurate 
phase described with v = (0, 0, 0) for both the Ni s and Ni c 
sites. The components of the vector correspond to the spin 
component on the Ni sites in the order given in Table fTTTl 
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TABLE VI: Irreducible representation of the group G v for the 
incommensurate magnetic structure with v = (9, 0,0). 



space group: 

{1, 2 X , m xyi m xz } , (A2) 

The group consists of 4 different classes and therefore has 
4 one-dimensional irreps. We determined the classes and 
the character table of this group and these are shown in 
Table IVT1 

The eigenvectors were calculated using the projector 
method using Eq. IjAljl and are given in Table IVIll 

As explained in the text, we replace the variables in the 
above table by symmetry adapted coordinates, so that we 
parametrize the eigenvectors as in Table IVIIll 

To fit to data one assumes a single irrep and then op- 
timizes with respect to the choice of the complex- valued 
amplitudes m°, m b s , etc. For the HTI phase one selects 
the representation which best fits the data for the opti- 
mized values of the amplitudes. Alternatively, one can 
introduce the symmetry adapted coordinates as in Eqs. 
© and accept the conclusion from Landau theory that 
these symmetry adapted coordinates all have the same 
complex phase. These symmetry adapted coordinates 
are constructed by the methods leading to Eqs. JSJ and 
are given in Table IVnTl 
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TABLE VII: Irreducible representations for the incommensu- 
rate phase associated with v = (q, 0, 0) for the Ni s and Ni c 
sites. The components of the vector correspond to the spin 
component on the Ni sites in the order given in Table UTTl 



through 36 

W) = f 2 ?) 2 N m ^$ R(Q) |/(Q)| 2 |F ± (Q)| 2 , 

\ £ ' V m o 

(Bl) 

where /(Q) is the magnetic form factor for Ni 2+ ions. 37 
<J> is the flux of incident neutrons, N m and V m o are 
the number and the volume of the magnetic unit cells, 
7 = -1.913 and r = 2.818 • KT 15 m. J(Q) is the to- 
tal integrated intensity of a Bragg reflection when the 
sample is rotated about the normal to the scattering 
plane and i?(Q) is a factor which takes into account 
the Q-dependent sensitivity of the spectrometer. For the 
powder experiment, R(Q)—1/ sin(6'(Q)). For the triple- 
axis experiment, i?(Q) was calculated using the Cooper- 
Nathans approximation. 38 F_l(Q) is the component of 
the magnetic structure factor perpendicular to the scat- 
tering wave vector and is defined as 

F ± (Q)=F(Q)-(F(Q)-Q)Q, (B2) 

where Q = Q/Q. 

The magnetic structure factor for a commensurate 
structure was calculated using 
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TABLE VIII: Symmetry adapted coordinates which trans- 
form according to the irreps for the incommensurate phase 
associated with v = (q, 0, 0) for the Ni s and Ni c sites. When 
more than one representation is active, it is a reasonable ap- 
proximation to allow all parameters of representation F n to 
have the same phase factor e"*" . 



APPENDIX B: MAGNETIC NEUTRON 
SCATTERING CROSS-SECTION 



F(Q) = ^m i exp(-iQ-d i ) (B3) 

i 

where are the positions of the Ni 2+ -ions and the sum 
is over the magnetic ions in the magnetic unit cell. The 
magnetic dipole moments are given by one of the 
irreps given in Table [V] For incommensurate structures 
the magnetic structure factor was calculated using 

F(Q) = ^^*y ex p(-iQ-d i ), (B4) 

i 

where the Fourier components of the structure, , are 
given by the basis vectors of the irreps of the little group 
given in Table I VIII The moment on site i is given by 
m i =~(\I> v + (\& v )*), so the factor \ in Eq. l|B4jl ensures 
that the scale factor for the basis vectors of the irreps 
can be expressed in Bohr magnetons. 

APPENDIX C: ABSOLUTE MAGNITUDE OF 
THE ORDERED MOMENT 

The magnitude of the ordered moments were deter- 
mined by comparing the nuclear and the AF Bragg peak 
intensities. The intensity of a nuclear Bragg peak is given 
as 

I(Q) = iV MC %^$ R(Q) \F N (Q)\ 2 , (CI) 

where N uc and Vo are the number and the volume of the 
unit cells. i*V(Q) is the nuclear structure factor given 
asS 



The integrated intensity of magnetic Bragg peaks is 
related to the structure factor of the magnetic ordering 



Fn(Q) =^6 i exp(-iQ-d i ), (C2) 

i 
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where the sum runs over all elements in the nuclear unit 
cell and bi is the bound coherent scattering length^ of 
atom i in the unit cell. Since multiple scattering and 
extinction corrections were important for the strongest 
reflections of the single crystal, their intensities were ig- 
nored for the normalization. The measured nuclear in- 



tensities gave the overall scale factor A cxp = N uc - V J $ 
m Eq The magnitude of the ordered magnetic mo- 

ment was determined using A Bxp and Eq. IjBlfl , allowing 
then the determination of the ordered moment for the 
various structures. 
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